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Abstract. Lie symmetry analysis is one of the powerful tools to analyze nonlinear 
ordinary differential equations. We review the effectiveness of this method in terms 
of various symmetries. We present the method of deriving Lie point symmetries, 
contact symmetries, hidden symmetries, nonlocal symmetries, A-symnretries, adjoint 
symmetries and telescopic vector fields of a second-order ordinary differential equation. 
We also illustrate the algorithm involved in each method by considering a nonlinear 
oscillator equation as an example. The connections between (i) symmetries and 
integrating factors and (ii) symmetries and integrals are also discussed and illustrated 
through the same example. The interconnections between some of the above 
symmetries, that is (i) Lie point symmetries and A-symmetries and (ii) exponential 
nonlocal symmetries and A-symmetries are also discussed. The order reduction 
procedure is invoked to derive the general solution of the second-order equation. 
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During the past four decades a great deal of interest has been shown in identifying and 
characterizing qualitative and quantitative properties of finite dimensional integrablc 
nonlinear dynamical systems B El- Several powerful mathematical methods have 
been developed to solve/integrate nonlinear ordinary differential equations (ODEs). 
Some widely used methods to solve nonlinear ODEs in the contemporary literature 
are (i) Painleve singularity structure analysis m H % (ii) Lie symmetry analysis 
0 m s u uni E], (iii) Darboux method [12] and (iv) Jacobi last multiplier method 
mmm- Among these the Lie group method advocated by Sophus Lie plays a vital 
role The method is essentially based on the invariance of differential 

equations under a continuous group of point transformations and such transformations 
usually have the form T = T(t,x,e), X = X(t,x,e), where (f,x) and (T, A") are the 
old and new independent and dependent variables respectively of the given ODE and 
e denotes the group parameter. The transformations depend only on the variables t 
and x and not on the derivatives, that is x. Transformations of this type are called the 
point symmetry group of a differential equation when this group of transformations leave 
the differential equation invariant mum- The Norwegian mathematician, Sophus 
Lie, founder of this method, developed an algorithm to determine the symmetry groups 
associated with a given differential equation in a systematic way. Once the symmetry 
group associated with the differential equation is explored, it can be used to analyze 
the differential equation in several ways. For example, the symmetry groups can be 
used (i) to derive new solutions from old ones UU, (h) to reduce the order of the given 
equation ME], (hi) to discover whether or not a differential equation can be linearized 
and to construct an explicit linearization when one exists mmm and (iv) to derive 
conserved quantities [6]. 

However, studies have also shown that certain nonlinear ODEs which are integrable 
by quadratures do not admit Lie point symmetries mu Eg. To understand the 
integrability of these nonlinear ODEs, through Lie symmetry analysis, attempts have 
been made to extend Lie’s theory of continuous group of point transformations in 
several directions. A few notable extensions which have been developed for this 
purpose are (i) contact symmetries [231 EH EE EH, (ii) hidden and nonlocal symmetries 
[36l EH EU E3 ESI E3 EH ESI EH [331 EH E5], ( ui ) A-symmetries [201 EH m SO], (iv) 
adjoint symmetries mu El HU and (v) telescopic vector holds mm- 

In the conventional Lie symmetry analysis the invariance of differential equations 
under one parameter Lie group of continuous transformations is investigated with point 
transformations alone. One may also consider the coefficient functions £ and t] (see 
Eq.(2.2)) in the infinitesimal transformations to be functions of x besides t and x. Such 
derivative included transformations are called contact transformations. In fact, Lie 
himself considered this extension JJ3j. The method of Ending contact symmetries for 
certain linear oscillators (harmonic and damped harmonic oscillators) were worked out 
by Schwarz and Cervero and Villarroel [21, Tl]. The integrability of a class of nonlinear 
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oscillators through dynamical symmetries approach was carried out by Lakshmanan and 
his collaborators, see for example Refs. [45] 46]. 

Investigations have also revealed that nonlinear ODEs do admit nonlocal 
symmetries. A symmetry is nonlocal if the infinitesimal transformations depend upon an 
integral. Initially some of these nonlocal symmetries were observed as hidden symmetries 
of ODEs in the following manner. Suppose an n th -order ODE is order reduced to (n—l) th 
order with the help of a Lie point symmetry. Now substituting the transformation (which 
was used to reduce the n t,l -order to (n — l) th order) in other symmetry vector fields of 
the n th -order equation one observes that these symmetry vector fields turns out to be 
symmetry vector fields of the order reduced ODE. In other words, all these vector fields 
satisfy the linearized equation of the order reduced equation. Upon analyzing these 
vector fields one may observe that some of them retain their point symmetry nature 
and the rest of them turn out to be nonlocal vector fields of the reduced ODE. Since 
these nonlocal symmetry vector fields cannot be identified through Lie point symmetry 
analysis they are coined as hidden symmetries. These nonlocal hidden symmetries were 
first observed by Olver and later they were largely investigated by Abraham-Shrauner 
and her collaborators [25], 2?, ]29] • 

Subsequently it has been shown that many of these nonlocal or hidden symmetries 
can be connected to A-symmetries. The A-symmetries concept was introduced by Muriel 
and Romero 05- These A-symmetries can be derived by a well-defined algorithm 
which include Lie point symmetries as a specific sub-case and have an associated order 
reduction procedure which is similar to the classical Lie method of reduction [38]. 
Although, A-symmetries are not Lie point symmetries, the unique prolongation of vector 
fields to the space of variables (t,x,x,x ,...) for which the Lie reduction method applies 
is always a A-prolongation, for some functions A (t,x,x,x,...). For more details on A- 
symmetries approach, one may refer the works of Muriel and Romero [39]. The method 
of finding A-symmetries for a second-order ODE has been discussed in depth by these 
authors and the advantage of finding such symmetries has also been demonstrated by 
them. The authors have also developed an algorithm to determine integrating factors 
and integrals from A-symmetries for second-order ODEs [38] . 

Very recently Pucci and Sacomandi have generalized A-symmetries by introducing 
telescopic vector fields. Telescopic vector holds are more general vector fields which 
encompose Lie point symmetries, contact symmetries and A-symmetries as their sub¬ 
cases. For more details about these generalized vector fields we refer to the works of 
Pucci and Sacomandi [ 42] . 

The connection between symmetries and integrating factors of higher order ODEs 
was investigated by several authors [47]. The literature is large and in this paper we 
discuss only one method, namely adjoint symmetry method which was developed by 
Bluman and Anco m EH- The main observation of Bluman and Anco was that the 
integrating factors are the solutions of adjoint equation of the linearized equation. In 
case the adjoint equation coincides with the linearized equation, then the underlying 
system is self-adjoint and in this case the symmetries become the integrating factors. 
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A main advantage of this method is that we can find the integrals straightaway by 
multiplying the integrating factors and integrating the resultant ODE [41] . 

The symmetry methods described above are all applicable to any order. Each 
method has its own merits and demerits. In this paper, we review the methods of 
finding symmetries (starting from Lie point symmetries to telescopic vector fields) of a 
differential equation and demonstrate how these symmetries are helpful in determining 
integrating factors and integrals of the ODEs and establish the integrability of them. 
We demonstrate all these symmetry methods for a second-order ODE. The extension 
of each one of this procedure to higher order ODEs is effectively an extension. We 
also illustrate each one of the methods with the same example. We consider the same 
example for all the methods so that the general reader can understand the advantages, 
disadvantages and complexity involved in each one of these methods. 

The example which we have chosen to illustrate the symmetry methods is the 


modified Emden equation (MEE) [1H1 H9j [501 ED ED EH [511 E5] 

o d 

x + 3xx + x' — 0, . = —. 

dt 


( 1 ) 


In the contemporary literature Eq. (JTJ) are also called second-order Riccati 
equation /Painleve-Ince equation. This equation has received attention from both 
mathematicians and physicists for a long time [561 E] • For example, Painleve had studied 
this equation with two arbitrary parameters, x + axx+/3x 3 = 0, and identified Eq. © as 
one of the four integrable cases of it [3], The differential equation ([I]) arises in a variety of 
mathematical problems such as univalued functions defined by second-order differential 
equations EH and the Riccati equation [48]. On the other hand, physicists have shown 
that this equation arises in different contexts. For example, it arises in the study of 
equilibrium configurations of a spherical gas cloud acting under the mutual attraction 
of its molecules and subjected to the laws of thermodynamics [501 SI]- Equation ((D) 
admits time independent nonstandard Lagrangian and Hamiltonian structures [53]. In 
the contemporary literature, this equation has been considered by several authors in 
different contexts. For example, Chandrasekar, Senthilvelan and Lakshmanan have 
studied the linearization and investigated the integrability of this equation through 
the extended Prclle-Singer procedure [53]. The Lie point symmetries of this equation 
were also derived by few authors in different contexts. For example, Mahomed and 
Leach have studied the invariance of this equation and shown that it admits sl(3, R) 
algebra and constructed a linearizing transformation from the Lie point symmetries. 
Pandey et al. have identified ([D) as one of the nonlinear ODEs that admits maximal 
Lie point symmetries when they carry out Lie symmetry analysis for the equation, 
x + f{x)x + g(x) = 0, where f(x) and g(x) are functions of x [58j. Nucci and her group 
have analyzed this equation in terms of Jacobi last multiplier [15]. A-symmetries and 
their associated integrating factors for this equation were investigated by Bhuvaneswari 
et.al. [59]. Noether symmetries of this equation were also studied in Ref. [60]. The 
nonlocal symmetries were also investigated in the Refs. mm 

The plan of the paper is as follows. In Sec. [21 we present Lie’s invariance analysis 
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to determine Lie point symmetries of Eq.dU). We also discuss few applications of Lie 
point symmetries. In Sec. [3l we describe the method of finding variational symmetries 
of this equation. In Sec. [2 we consider a more generalized transformation and present 
the method of finding contact symmetries of the given differential equation. In Sec. [5J 
we discuss the methods that connect symmetries and integrating factors. We consider 
two different methods, namely, (i) adjoint symmetries method and (ii) A-symmetries 
approach. In Sec. [6] we introduce the notion of hidden symmetries and enlist some 
hidden symmetries of MEE which are not obtained through Lie symmetry analysis. 
In Sec. [?1 we introduce nonlocal symmetries and investigate the connection between 
nonlocal symmetries and A-symmetries. In Sec. [HI we consider a more general vector 
held, namely telescopic vector held and derive these generalized vector helds for the 
MEE. Finally, we give our conclusions in Sec. [9] 

2. Lie point symmetries 

Let us consider a second-order ODE 

x = cf)(t, x, x), (2) 

where overdot denotes differentiation with respect to L t\ The invariance of Eq. (J2]l under 
an one parameter Lie group of inhnitesimal point transformations [6lj 158], 

T = t + e£(t,x), X — x + sr)(t, x), e«l, (3) 

where £(t,x) and rj(t,x ) are arbitrary functions of their arguments and e is a small 
parameter, is given by 

+ (Vt + x(r} x - &) - - (vtt + (27 ] tx - £ tt )x + (r) xx - 2£ tx )x 2 

- txxX 3 + (jh - 2ft - 3x£ x )x) = 0. (4) 

Substituting the known expression 0 in (J4J) and equating the coefficients of various 
powers of x, to zero one obtains a set of linear partial differential equations for the 
unknown functions f and r/. Solving them consistently we can get the Lie point 
symmetries (f and rj) associated with the given ODE. The associated vector held is 
given by V = f f + rj£. 

One may also introduce a characteristics Q — r) — xC, and rewrite the invariance 
condition ([4]) in terms of a single variable Q as [6j 

M = 0. (B) 

dt 2 dt y ^ y ' 

Solving Eq.(j5]) one can get Q. From Q one can recover the functions f and rj. The 

invariants associated with the vector held V can be found by solving the following 

characteristic equation m 

dt dx dx 

f V 1? (1> ' 


( 6 ) 
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Here r/ ]) represents the first prolongation which is given by ry + x{r) x — £ t ) — x 2 £ x . 
Integrating the characteristic equation dSJ) we obtain two invariants, namely u(t,x ) and 
v(t,x,x). The derivative of these two invariants, 
dv v t + xv x + 4>v± 

w = T =--> ( 7 ) 

du Uf + xu x 

defines a second-order differential invariant. Integrating the above equation ([JJ) we can 
get the solution for the given equation (J2j). 


2.1. Example: modified Emden equation 

Eq.flT]) is invariant under the infinitesimal transformation (J2J) provided it should satisfy 
the Eq.(j4j). Substituting the expression = —3xx — x 3 in (@|, we get 

r]{—3x - 3x 2 ) + (rj t + xrj x - x£ t - x 2 £ x )(-3x) - {r) tt + x{2r) tx - £«) 

+ i 2 (Vxx ~ 2 £tx) - ZxxX 3 + {rjx ~ 26 - 3x£ x )(-3xx - x 3 )) = 0. (8) 

Equating the coefficients of various powers of x to zero and solving the resultant set of 
partial differential equations for £ and rj, we obtain 

°2 + bi ^ 2 


2 

3 


Cl + do 3 di 4 

2 4 


£ = x(a 2 + afi - 

- y ^ 3 + (ci + 7}d 2 ^t 2 + feit + fe 2 , 

3 f , . , d\ 4 [Ci + d 2 \ 3 C 2 + fel 2 

V = ~ x [afi + a 2 + —t - —-—Jt -- —t 


+ x 2 ( dfi 3 
3 


Ci + d 2 


' ^t 2 — (c 2 + fei)t + Gq^ 


+ x^ — — dit 2 + Cit + c^j + dfi + d 2 , 


(9) 


where ai,bi,Ci and di, i = 1,2, are real arbitrary constants. The associated Lie vector 
fields are given by 


Ei = 


d 


T j. . ^ xt\ d 9 / xt\ d 

V2 = tn ~T)m +xt \~ 1 + Y)di' 


dt ’ 

E = x-^- — x 3 -^-, V 4 = xt^— + x 2 ( 1 


dt 
xt 2 d 


dx dt 

„„ „ . x 2 t 2 \ 9 

1/5= --s + d 1 ^ :c(+ —)&■ 


. 9 


E> = t^l 
E = ^d 


x 2 t 2 \ 9 


xt\ 9 / 3 

— ITT+Xtll-Xt + , „ , 

2/at v 2 2 yax’ 


xt\ d 

J )a + |1 


3 2 2 x 3 t 3 \ d 

-xt 


t . t 3 / xt\ d / 3 29 

Vs = -fi 1 -Y)m +t { 1 -2 xt + xt 


x 3 t 3 \ d 


A Jdx ^ 

One can explore the algebra associated with the Lie group of infinitesimal 
transformations (j3J) by analyzing the commutation relations between the vector fields. 
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Since the example under consideration admits maximal Lie point symmetries (eight) 
the underlying Lie algebra turns out to be sl( 3, R) which can be unambiguously verified 
from the vector fields m [18]. In the following, we present few applications of Lie 
vector fields (fTOl) . 


2.2. Applications of Lie point symmetries 

2.2.1. General solution The first and foremost application of Lie point symmetries is 
to explore the solution of the given equation through order reduction procedure. The 
order reduction procedure is carried out by constructing the invariants associated with 
the vector fields [ 6 : 9J. In the following, we illustrate the order reduction procedure by 
choosing the vector field V 3 . For the remaining vector fields one can proceed and obtain 
the solution in the same manner. 

Substituting the expressions £, and rj^ in the characteristic equation j = 
Jy, we get 


dt dx 


dx 


x 


-x c 


-{3xx 2 + x 2 )' 

Integrating the characteristic equation dll]) we find the invariants u and v are of the 
form 


u = t 


v = 


x(x + x 2 ) 


X 


X 


( 12 ) 


The second-order invariant can be found from the relation w = ^ (see Eq. flTjl ). 
Substituting Eq. (1T2T) and their derivatives in (J7J) and simplifying the resultant equation 
we arrive at 


2 4 

X X 


^ = X 2 (1 + 2-+ ;_) = !> 2 . 


du 


x x* 


(13) 


Integrating equation m we find v = where I\ is an integration constant. 

Substituting the expressions u and v in this solution and rewriting the resultant equation 
for x, we end up with 

h+t ^ J 

Integrating Eq. (11411 we obtain the general solution of the MEE equation in the following 
form 

2(/i + 1 ) 


x(t) = 


(15) 


I 2 + 2I\t + t 2 

where 1 2 is the integration constant. The solution exactly coincides with the one found 
by other methods 
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2.2.2. Linearization One can also identify a linearizing transformation from the Lie 
point symmetries if the given equation admits sl(3,R) algebra [18]. The method of 
finding linearizing transformation for the modified Emden equation was discussed in 
detail by Mahomed and Leach [61] and later by others [62] [53]. The underlying idea 
is the following. One has to choose two vector fields appropriately in such a way that 
they should constitute a two-dimensional algebra in the real plane [33] and transform 
them into the canonical form and tjr, where t and x are the new independent and 
dependent variables. For the MEE one can generate this subalgebra by considering the 
vector fields Vg and Vg = V? — 2Vq. To transform them into the canonical forms one 
should introduce the transformations [[T5] 

- tX ^ 1 (16) 


t = 


X = 


t(tx — 2 )’ 2t(tx — 2) J 

where t and x are the new independent and dependent variables respectively. In these 
new variables, (t,x), Eq.(pQ) becomes the free particle equation (£§ = 0. From the 
solution of the free particle solution one can derive the general solution of MEE ([!]) with 
the help of (TT 6 T 1 . The solution coincides with the one given in Eq. (fT5jh 


2.2.3. Lagrangian from Lie point symmetries Another interesting application of Lie 
point symmetries is that one can explore Lagrangians for the given second-order 
differential equation through the Jacobi last multiplier [15]. The inverse of the 
determinant A m, 


A = 


1 

6 


X 

X 


( 1 ) 

Vi 

Vi 


(i) 


m 


M=-, 


(17) 


where (£ 1 , ?/i) and (£ 2 , V 2 ) are two sets of Lie point symmetries of the second-order ODE 
([2]) and 77] 11 and are their corresponding Erst prolongations, gives the Jacobi last 
multiplier for the given equation. The determinant establishes the connection between 
the multiplier and Lie point symmetries. The Jacobi last multiplier M is related to the 
Lagrangian L through the relation [14] 

d 2 L . , 

(18) 


M = 


dx 2 ' 


Once the multiplier is known, the Lagrangian L can be derived by integrating the 
expression w two times with respect to x. 

To obtain the Jacobi last multiplier M for the MEE, we choose the vector fields V 3 
and V\. Evaluating the determinant (1171) with these two vector fields, we hnd 

—3xx — x 3 


so that 


A = 


M = 


1 

x 

1 


X 
—X' 

0 


3 —x(x + 3x 2 ) 

0 


= —(x 2 +1) 3 


(X“ 


+ ±y 


(19) 


( 20 ) 
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We can obtain the Lagrangian by integrating the expression (I2U1) twice with respect to 
x. Doing so, we find 

L= ~ 2(x + x 2 ) + f 1 ( t ’ x )x + f2(t,x), ( 21 ) 

where fi and / 2 are two arbitrary functions. The Lagrangian leads to the equation of 
motion ([T]) with One can also find more Lagrangains for Eq.((2]) by considering 

other vector fields given in (fTUj) . 

3. Noether symmetries 

In the previous section, we have discussed the invariance of the equation of motion under 
the one parameter Lie group of infinitesimal transformations (|3]). If the given second- 
order equation has a variational structure then one can also determine the symmetries 
which leave the action integral invariant. Such symmetries are called variational 
symmetries. Variational symmetries are important since they provide conservation laws 
via Noether’s theorem [63]. In the following, we recall the method of finding variational 
symmetries mm- 

Let us consider a second-order dynamical system described by a Lagrangian 
L(t,x,x ) and the action integral associated with the Lagrangian be 

S — j L(t,x,x)dt. (22) 

Noether’s theorem states that whenever the action integral is invariant under the 
one-parameter group of infinitesimal transformations (]3]) then the solution of Euler’s 
equation admits the conserved quantity [6l, 65] 

81 

/= ( ^-77)—-£L + /, (23) 

where / is a function of t and x. The functions £, ij and / can be determined from the 
equation 

= € + ’’§§ + «-*>ir (24) 

where over dot denotes differentiation with respect to time and 

E{L} = -iL + f. (25) 

Equation (124|) can be derived by differentiating the equation il23l) and simplifying the 
expression in the resultant equation. Solving equation (|24l) one can obtain explicit 
expressions for the functions £, 77 and /. Substituting these expressions back in (1231) 
one can get the associated integral of motion. 

3.1. Example: modified Emden equation 

In this sub-section, we illustrate the method of finding Noether symmetries and their 
associated conserved quantities for the MEE (JT[) which has a nonstandard Lagrangian 
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of the form (see Eg. (1211) ) 

L= 3{x + x 2 Y (26) 

where we have chosen the arbitrary functions f\ and f 2 to be zero for simplicity. 
Substituting the Lagrangian (1261) and its derivatives in (1241) . we get 

2x 


n\ - 


3(x + x 2 ) 2 

= -«.+*«.) ( 3 n 


^ + ft + Xf x - 


Of • 1 9\9 / 1 JT- ' ' ) 

,3{x + x 2 ) 2 / 

Equating the coefficients of various powers of x to zero and solving the resultant 
equations, we find 

77 = 12 D — 6 (C + 4 Dt)x + 3 (B + 3 t(C + 2 Dt))x 2 
+ t{2B + 3Ct + Wt 2 ))x 3 , 

£=E- 3t{C + 2 Dt) + ^(4 + t( 2 B + 3 Ct + 4Dt 2 ))x, 
f = At+ Bt 2 + Ct 3 + Dt 4 , (28) 

where A, B , C, D and E are real arbitrary constants. The associated vector fields are 


= 4 - * 4 ' A ' 2 = 4 + 1 


— tx 


a 


J dx' 


3 t 2 x\ d 


2 J dt 

_ /t 3 x t 2 \d 


3 t 2 x 3 \ d 


X 3 = [t -— — + [2x - 3tx~ H-— — 


2 tx + 


2 

3 t 2 x 2 


J dx ’ 
t 3 x 3 \ 9 


d 


2 J dx 1 dt' 


(29) 


The Noether’s symmetries are sub-set of Lie point symmetries. In the above, while 
the vector fields X l 7 X 2 and X 5 exactly match with the Lie vector fields V 3 , V4 and V\ (see 
Eq. (PHD ), the remaining two Noether vector holds X 3 and X 4 can be expressed as linear 
combinations of other Lie point symmetries, that is A " 3 = V 2 + 2 V 5 and X 4 — V? — 2Vq. 

Substituting each one of the vector fields separately into (1231) we obtain the 
associated integrals of motions. They turned out to be 

(—x + tx 2 + tx ) 2 


h=t- 
h = 
h = 
h = 


x 


h = 


X 2 + X ( X 2 + x ) 2 

—9 t 2 x 3 + 3t 3 a : 4 — 3x(2 + 3 t 2 x) + 2tx 2 (6 + 3 t 2 x) + 9tx(6 + 3 t 2 x) 


(. x 2 + x) 2 


3(2 — 2 tx + t 2 x 2 + t 2 x) 


(. x 2 + x) 


2x + x 2 


dh 


= 0 , % = 1 , 2 , 3 , 4 , 5 . 


(30) 


3{x 2 + x) 2 ' dt 

One can easily verify that out of the five integrals two of them are independent 
and the remaining three can be expressed in terms of these two integrals, that is 
I 2 = If, I 3 = / 1/4 and J 5 = I(I 4 — 3 If). We can construct a general solution of 
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dH) with the help of the two independent integrals A and J 4 . The underlying solution 
coincides with (fl5l) after rescaling. 


4. Contact symmetries 

In the previous two cases, Lie point symmetries and Noether symmetries, we have 
considered the functions £ and 7 ] to be functions of t and x only. One may relax this 
condition by allowing the functions £ and 7 ] to depend on x besides t and x. This 
generalization was considered Sophus Lie himself [43] and later by several authors 
[2U E2l sa 06]. This velocity dependent infinitesimal transformations are called 
contact transformations and the functions £ and t] are called contact symmetries. The 
contact symmetries for the harmonic and damped harmonic oscillators were worked out 
explicitly by Schwarz and Cervero and Villarroel respectively [2T1 22 \. Several nonlinear 
second-order ODEs do not admit Lie point symmetries but they were proved to be 
integrable by other methods. To demonstrate the integrability of these nonlinear ODEs 
in Lie’s sense one should consider velocity dependent transformations. In the following, 
we give a brief account of the theory and illustrate the underlying ideas by considering 
MEE as an example. 


4-1. Method of Lie 


Let a one-parameter group of contact transformation be given by 

T — t + e £(£, x, x), X — x + £T]{t, x, x), X — x + £ //^(f, x, x) e-C 1. (31) 

The functions £ and t] determine an infinitesimal contact transformation if it is possible 
to write them in the form 


f(t,x,x) = 


dW 


7 ]{t , X, x) = W — X 


dW m dW . dW 
- T — ——— x - 


(32) 


dx dx dt dx 

where the characteristic function W ( t , x, cc) is an arbitrary function of its arguments. 
If W is linear in x the corresponding contact transformation is an extended point 
transformation and it holds that W(t,x, x) = 7 ](t,x ) — xt;(t,x). A second-order 
differential equation ([2]) is said to be invariant under the contact transformation (|3TT) 
if + Vfa + ? 7 ^§§ — — 0 on the manifold x — <f(t,x,x) = 0 in the space of the 

variables t,x,x and x [21, 22]. Here and r/ 2 ^ are the first and second prolongations 
with ?/b = fj — xf and r — xf. The invariance condition provides the following 
linear partial differential equation for the characteristic function W: 


dW df 
dx dt 


+ x- 


dW 

dx 



'dW . dW 

~W + X ~d. 7 



2 d 2 W nl d 2 W 
+ 20 - 


dx 2 


dtdx 


0 , . d 2 W 
+ 2 cj)x n n . + 


dW d 2 W d 2 W 

+ + 2X- 


2 d 2 W\ 

+ x T^)-°- 


dxdx r dx dt 2 1 dtdx 1 ~ ' 

one can get the characteristics W. From W one can recover the 
contact symmetries £ and 77 . One can also recover the necessary independent integrals 
from the characteristic function (see for example, Ref.[ 


Integrating Eq. 
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To determine the contact symmetries of MEE, we have to determine the characteristic 
function (l33]h by solving the first-order partial differential equation 


(3i + 3x 



3\2 


(3xx + x 6 ) 


d 2 W 


dx 2 


— 2(3xx + x s ) 


dtdx 


dxdx 

d 2 W n . d 2 W 
2 x- 


(3xx + x d ) 


dW 


dx 


+ x 2 


d 2 W 


= 0 . 


dt 2 ' ~~ dtdx ' ~ dx 2 ) 

One may hnd the general solution of the above linear partial differential equation 
by employing the well known methods for solving linear partial differential equations. In 
general W depends upon arbitrary functions and the contact Lie group has an infinite 
number of parameters. 

For the sake of illustration, in the following, we present two particular solutions of 

Eq.flMD: 

X (x 2 + x) 


W 1 = 


x 2 (l 


tx) - fi, 


X 2 + X 

The infinitesimal vector fields read 

d x 2 {l — tx) d 


W 2 = 


xx 


\Jx 2 + 2 x \Jx 2 + 2 x 


— t — + 
dt 


x + x 2 dx ’ 


— 


x 


d {x + x 2 ) d 


x- 


y/ 2 x + x 2 dt y/ 2 x + x 2 dx 


(35) 


(36) 


As one can see from (1361) the infinitesimals £ and 77 depend on the velocity terms 
also. One can derive the general solution from each one of the contact symmetries 
by solving the characteristic equation associated with the vector field. We demonstrate 
this procedure in Sec. 14.31 


4-3. Method of Gladwin Pradeep et.al 

In a recent paper Gladwin Pradeep et.al proposed a new method of finding contact 
symmetries for a class of equations [23]. Their method involves two steps. In the first 
step, one has to hnd a linearizing contact transformation. In the second step, with 
the help of contact transformation, one proceeds to construct contact symmetries for 
the given equation. Once the contact symmetries are determined the order reduction 
procedure can be employed to derive the general solution of the given differential 
equation. In the following, we recall this procedure with MEE as an example. 

Step 1: Linearizing contact transformation 

The MEE (JT]) can be linearized to the free particle equation = 0 by the contact 
transformation (for more details one may refer, Ref. 


x = 


2 uu 
1 + u 2 ’ 


X = 


2w 2 (l 


u 


(1 + u 2 ) 2 


(37) 
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One may also invert the above transformation (1HTD and obtain 


x 


u = 


y/ 2 x + x 2 


u = 


(x + x 2 ) 
\j 2 x + x 2 


(38) 


Step 2: Contact symmetries 

Let us designate the symmetry vector field and its first prolongation of the nonlinear 


ODE (pQ) respectively be of the form D = \4- + and 


9 


d 


dt 

d 




(39) 


dt dx dx' 

where A and /i are the infinitesimals associated with the variables t and x, respectively. 

Let the symmetry vector field associated with the linear ODE, = 0, be 
A = and its first extension be A 1 = + 17 -^ + (?) — m£)JU Using the 

contact transformation ()38li we can deduce the following differential identities, that is 

d_ 

du du dx 1 du dx W,J ” LVA Qu 

these relations, we find 


9x 9 1 9± 9 and 737 : = Rewriting the first prolongation A 1 using 


AI = 4 + (”i + «-^il& 


d (dx • dx 

+ I v— + [V ~ «?) 


Now comparing the vector fields (l40]i and 


A = & 


dx 


M= „_ + (li _ o- = 


9m 

, we find 
9x \/ 2 x + x 2 


_ d_ 

du ) d x 


{rjx + rjx) — £x. 


(40) 


(41) 


X + X 2 

The functions £ and rj are the symmetries of the free particle equation. They are given 
by m 


Ai = 
A fi = t- 


Ao — 


d 


du ’ 


. d 
As — 

du 
d 


A d = u 


d 


du ’ 


A 5 = w 
9 


9 


90 


(42) 


d_ 
dt J 

d 2 9 9 9 

'S’ A7 = t a +i “S? As = i “a + “'S? 

Substituting the above symmetry generators Aj’s, i = 2,..., 8 , in Eq. (lTT]h we can 
determine the function /i. Substituting A = £ and // in the vector field D = £f| + 
we arrive at the following contact symmetry generators of Eq. ([I]), 


_ 9 V 2 x + x 2 . 9 

Dl — W7) D 2 — -—-^ "3 — 

9r a; + x z 9x 


\/ 2 x + x 2 . . . 9 

• , 2 ( tx + x )a~» 

x + x z 9x 


D 4 = 


x( 2 x + x 2 ) 9 
(x + x 2 ) 9x’ 


Dr — 


X 


9 


(x + X 


2 ^ 5 


x- 


9 9 2 „ 

a ^ t dt ~ z ai' Slj = t at + 


\/2x + x 2 9f o/ 2 x + x 2 9x ’ 
9 x 2 (l — tx) 9 


Da — 


tx 


9 


+ 


x + x 2 dx ’ 
x 2 y / ( 2 x + x 2 ) fx(x + x 2 ) \ 9 


y/ (2x + X 2 ) ^ ' X + X 2 y/ (2x + X 2 ) 


9x 


(43) 


One can unambiguously verify that all the symmetry generators are solutions of the 
invariance condition. 
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4-3.1. General solution To derive the general solution of the given equation one has 
to integrate the Lagrange’s system associated with the contact symmetries given above. 
We demonstrate this procedure by considering the vector held ^4 given in Eq. (l43lh For 
the other vector fields one may follow the same procedure. 

The characteristic equation associated with the vector held hR turns out to be 
dt (x + x 2 )dx (x + x 2 )dx 

0 x(2x + x 2 ) x A + x 2 x — 2x 2 ' 


Integrating Eq. (l44|i we find the invariants u and v are of the form u = t and v — - + x. 
In terms of these variables one can reduce the order of the equation (JTJ) . The reduced 
equation turns out to be the Riccati equation, ^ = — v 2 , whose general solution is given 
by v = where I\ is the integration constant. Substituting the expressions u and v 


I\-\-U ? 

in the above solution and rearranging the resultant expression for x, we find 
x 


x — 


I\ +1 


+ x = 0 . 


(45) 


Integrating (l45|i one can obtain the general solution of ([T]) . The general solution coincides 
with (fUjl) after appropriate rescaling. 


5. Symmetries and integrating factors 

In SecEl we have discussed only few applications of Lie point symmetries. One can 
also determine integrating factors from Lie point symmetries. I 11 fact, Lie himself 
found the equivalence between integrating factors and Lie point symmetries for the 
hrst-order ODEs. For second-order ODEs the equivalence has been established only 
recently [201]. The reason is that unlike the hrst-order ODEs (which admit infinite 
number of symmetries) several second-order ODEs do not admit Lie point symmetries 
although they are integrable by quadratures. Subsequently attempts have been made to 
generalize the classical Lie method to yield nontrivial symmetries and integrating factors. 
Two such generalizations which come out in this direction are (i) adjoint symmetries 
method pJDjHTj and (ii) A-symmetries approach [381 [39] . The adjoint symmetry method 
was developed by Bluman and Anco [9] and the A-symmetry approach was floated by 
Muriel and Romero. The applicability of both the methods have been demonstrated 
for the equations which lack Lie point symmetries. In both the methods one can find 
symmetries, integrating factors and integrals associated with the given equation in an 
algorithmic way. 

In the following, we recall briefly these two powerful methods and demonstrate the 
underlying ideas by considering MEE as an example. 

5.1. Method of Bluman and Anco 

In general, an integrating factor is a function, multiplying the ODE, which yields a first 
integral. If the given ODE is self-adjoint, then its integrating factors are necessarily 
solutions of its linearized system ([5]) [10]. Such solutions are the symmetries of the 
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given ODE. If a given ODE is not self-adjoint, then its integrating factors are necessarily 
solutions of the adjoint system of its linearized system. Such solutions are known as 
adjoint symmetries of the given ODE [9]. 

Let us consider second-order ODE ([2]). The linearized ODE for Eq. Q is given in 
(J5]). The adjoint ODE of the linearized equation is found to be 

. d 2 w d . , . i . 

L [x]w = w) - (j) x w = 0. (46) 

The solutions w = A (t,x,x) of the above equation (1461) holding for any x(t) satisfying 
the given equation (J2J) are the adjoint symmetries of (]2j) [9]. 

The adjoint symmetry of the Eq. (l2jl becomes an integrating factor of (J2J) if and only 
if A (t,x,x) satisfies the adjoint invariance condition m 

L* [x]A(t, x, x) = -A„(i - <p) + 4(A*(x -<p)). (47) 

Now comparing the Eqs. (1461) and (1471) and collecting the powers of x and constant 
terms, we find 

A t± + A x ±x + 2A X + A(j) xx + 2cf) x A x + (j> A xx = 0, (48 a) 

A t t T 2,Ai x x T A xx x -(- A (pt x T A (p xx x -)- (fi x A^ (p x A x x 
- A(p x - A x cj) + (f)A tx + 4>A xx x + A ± 4> t + A ± 4> x x = 0. (48 b) 

The solutions of (148 61) are called the adjoint symmetries. If these solutions also satisfy 
Eq. (148 aj) then they become integrating factors for the given second-order ODE (J2j). 

The main advantage of this method is that if the given equation is of odd order or 
does not have variational structure, one can use this method and obtain the integrals 
in an algorithmic way. 

5.1.1. Exainple: modified Emden equation For the MEE equation, the linearized 
equation is given by 

+ 3 x “7 T + ( 3;i; + 3x 2 )Q = 0. (49) 

dE dt 

The adjoint equation for the above linearized equation turns out that 

d 

—r + —(—3 xw) + (3x + 3 x 2 )w = 0. (50) 

dt z dt 

Since Eqs. (149]) and (150]) do not coincide, the function Q is not an integrating factor 
for the MEE. In this case the integrating factors can be determined from the adjoint 
symmetry condition (]50]) . The adjoint symmetry determining equation (1501) for the 
present example reads 

A tt + 2A tx x + A xx x 2 — 3Ax — 3xA t — 3xA x x + (3x + 3x 2 )A 
+ (3xx + x 3 )A x — (3xx + x 3 )A xt — ( 3xx + x 3 )xA xx — {3x + 3x 2 )A±x = 0. (51) 

Two particular solutions of (T5Tj) are given by 

,\ = x A = t{—2 + tx) 

1 (i + x 2 )’ 2 2 (tx — x + tx 2 ) 2 


( 52 ) 
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These two adjoint symmetries, A] and A 2 , also satisfy the Eg. (148 np . So they become 
integrating factors for the Eq.([T]). Multiplying the given equation by each one of these 
integrating factors and rewriting the resultant expression as a perfect derivative and 
integrating them we obtain two integrals I\ and J 2 which are of the form, 

(xt — x + tx 2 ) (2 + xt 2 — 2 tx + t 2 x 2 ) 


h = 


X + X 2 


h = 


2 {xt — x + tx 2 ) 

From these integrals, I\ and J 2 , the general solution can be derived, 
expression coincides with Eq. dlSp after rescaling. 


(53) 

The resultant 


5.2. Method of Muriel and Romero 

As we mentioned at few places earlier, many second-order nonlinear dynamical systems 
often lack Lie point symmetries but proved to be integrable by other methods. 
To overcome this problem, efforts have been made to generalize the classical Lie 
algorithm and obtain integrals and general solution of these nonlinear ODEs. One 
such generalization is the A-symmetries approach. This approach was developed by 
Muriel and Romero [2D]. These symmetries are neither Lie point nor Lie-Backlund 
symmetries and are called A-symmetries since they are vector fields that depend upon a 
function A. If we choose this arbitrary function as null we obtain the classical Lie point 
symmetries. The method of finding A-symmetries for a second-order ODE has been 
discussed in depth by Muriel and Romero and the advantage of finding such symmetries 
has also been demonstrated by them [38]. The authors have also developed an algorithm 
to determine integrating factors and integrals from A-symmetries for the second-order 
ODEs [39]. 

Consider a second-order ODE ([2]). Let V = £(f, x)-^ + r/(t,x be a A-symmetry 
of the given ODE for some function A = A (t,x,x). The invariance of the given ODE 
under A-symmetry vector field is given by |20| 

£4>t + V<t>x + - ry (2)l = 0, (54) 

where r/L-l 1 )] an d ^IM 2 )] are the first and second A- prolongations. They are given by 
, ; ^[A,(i)] _ _|_ X)rj^°^(t, x) — ( D t + A )(f(t, x))x = r+ A (r/ — ££), 

, ; ?[A,(2)] = _|_ a )'ij [x ’ {1)] (t,x) - (D t + A )(f(t,x))x = r/ (2) + /(A), (55) 

where /(A) is given by D[\](r) — fx) + 2A(r/d) — ££) + A 2 (?/ — fx) and r= rj(t, x), and 
D is the total differential operator (D = + xf^ + $-§£). In the above prolongation 

formula if we put A = 0, we end up at standard Lie prolongation expressions. Solving 
the invariance condition (154|) we can determine the functions £, r) and A for the given 
equation. We note here that three unknowns £, r/ and A have to be determined from 
the invariance condition (l54lh The procedure is as follows. 

Let us suppose that the second-order equation ([2]) has Lie point symmetries. In 
that case, the A-function can be determined in a more simple way without solving the 
invariance condition (l54]h 
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If V is a Lie point symmetry of ([2]) and Q = rj — x£ is its characteristics then v = 
is a A-symmetry of (J2]) for A = The A-symmetry satisfies the invariance condition 

m 

(f> x + A cf)± = D[X] + A 2 . (56) 

Once the A-symmetry is determined, we can obtain the hrst integrals in two different 
ways. In the hrst way, we can calculate the integral directly from the A-symmetry using 
the four step algorithm given below. In the second way, we can find the integrating 
factor fi from A-symmetry. With the help of integrating factors and A-symmetries we 
can obtain the hrst integrals by integrating the system of Eqs. (l58jh In the following, we 
enumerate both the procedures. 

(A) Method of finding the first integral directly from X-symmetry ISRl 
The method of finding integral directly from A-symmetry is as follows: 

(i) Find a hrst integral w(t, x, x ) of that is a particular solution of the equation 

w x + A w± = 0 , where subscripts denote partial derivative with respect to that 
variable and is the hrst-order A-prolongation of the vector held v. 

(ii) Evaluate D[w\ and express D[w\ in terms of (t,w) as D[w\ = F(t,w). 

(iii) Find a hrst integral G of d t + F(t,w)d w . 

(iv) Evaluate I(t,x,x) = G(t,w(t,x,x)). 

(B) Method of finding integrating factors from X Vfflf 

If V is a Lie point symmetry of (0|) and Q = g — xf is its characteristics, then 
v = d x is a X-symmetry of (dj) for X = D[Q\/Q and any solution of the first-order linear 
system 

m + (0i - = 0, = o, (57) 

is an integrating factor of (G|). Here D represents the total derivative operator and it is 
given by + 

Solving the system of equations (1571) one can get //. Once the integrating factor /i 
is known then a hrst integral / such that I± — /i can be found by solving the system of 
equations 

I t = g(Xx - 0), I x = -Xg, h = H- ( 58 ) 

From the hrst integrals, we can write the general solution of the given equation. 

5.2.1. Example: modified Emden equation Since the second-order ODE under 
investigation admits Lie point symmetries one can derive the A-symmetries directly 
from Lie point symmetries (59] through the relation A = 

To start with, we consider the vector held V 3 . In this case, we have £ = x and 
7 ] = —x 3 , the Q function turns out to be Q = 7] — £x = —x(x + x 2 ). Using the relation 
A = we can hx A 3 = | — x. In a similar way one can hx the A-symmetries for the 
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remaining vector fields. The resultant expressions are given in the following Table. One 
can verify that the functions £, r) and A satisfy the invariance condition (l54]h 


Vector 

Q 

A-symmetries 

Ui 

—x 

-(3 x + ^) 

d 2 

\t {-2 + tx)(x + x 2 ) 

2-xt z -itx+t z x' z 

t(2—tx) 

u 3 

—x(x + x 2 ) 

*-x 

X 

u 4 

x(—xt + x — tx 2 ) 

*-x 

X 

U 5 

|(2 + xt 2 — 2 tx + t 2 x 2 ) 

- — X 

X 

V 6 

|(—2 + tx)(xt — x + tx 2 ) 

f 2—xi‘‘—4tx+i' s x‘ 1 \ 

\ t(2—tx) > 

v 7 

|(2 — 3 t 2 x 2 + t 3 x 3 — 3 xt 2 + xt 3 x 

t(—x z t z +x(<o—§tx)-\-x z (6— §tx-\-t z x z ') 

2-3t 2 x 2 +t 3 x 3 +xt 2 (-3+tx) 

^8 

^(—2 + tx)( 2 + xt 2 — 2 tx + t 2 x 2 

2—xt 2 —4:tx+t 2 x 2 

t(2—tx) 


Table 1 . A- symmetries for eight the vector fields admitted by Eq. m 


Since we are dealing with a second-order ODE two different A functions are 
sufficient to generate two independent first integrals and hence the general solution. 
In the following, we consider the vector fields V 3 and Vg and their A-symmetries and 
demonstrate the method of finding their associated integrals. 


(i) First integrals directly from A- symmetry 


Substituting A 3 = ~ — x in the equation w x + Xw± = 0, one gets w x + — x^j w± = 0. 
This first-order PDE admits an integral w(t,x,x ) of the form w(t,x,x) = - + x (first 
step). The total differential of this function can be expressed in terms of w itself, that is 
D[w] = —w 2 = F(t,w ) (second step). Next, one has to find an integral associated with 
the first-order partial differential equation ^ — w 2 ^f = 0 . A particular solution of this 
first-order partial differential equation can be given as G = t — A (third step). Finally, 
one has to express G(t, w) in terms of (t, x, x). Doing so, we find that the integral turns 
out to be (fourth step) 


h=t- 


x 


d Ji = 0. 

dt 


(59) 


1 9 * 

x + x- 

Next we consider the function Ag. Following the steps given above, we find the 
integral associated with A 6 which turns out to be 


h = 


xt — X + tx 2 


2 + xt 2 — 2 tx + t 2 x 2 


(60) 


with ^ = 0. Using these two integrals, ii and J 2 , one can construct the general solution 
of Eq. ffTj) . The resultant expression coincides with the earlier expression (see Eq. (TT5T) ) 
after rescaling. 


(ii) Integrating factors from \-symmetries 
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In the following, we discuss the second route of obtaining the integral from A- 
symmetry. We solve the system of equations (j57j) in the following way. We first consider 
the second equation in (j571) and obtain a solution for p. We then check whether the 
obtained expression satisfies the first equation or not. If it satisfies then we treat it as 
a compatible solution. 

We again consider the Lie point symmetries V 3 and Vg and discuss the method 
of deriving integrating factors for these functions. To determine the integrating factor 
associated with A 3 we first solve the second equation in (1571) . that is + x)/j,x+^ti = 
0. A particular solution is //1 = — • This solution also satisfies the first equation in 


(1571) . To determine the integral we substitute pi and A 3 in 
set of equations for the unknown /, namely 

T — I t _ (* ~ ^ T . 

J-t -*-•) ± X ' - n t 1 nr. 


X 


and obtain the following 


(61) 


(x + x 2 ) 2 ’ (x + X 2 ) 2 

The integral which comes out by integrating the system of equations f| 6 ip . that is 
I\ — t — x 2 +± , coincides exactly with the one found earlier. 

Let us now consider the function X 6 . Substituting the function Ag in equation 
we get 

, 2 — xt 2 — 4 tx + t 2 x 2 


Eq. 


h* + ( t{2 _ tx) »)& - 0. 

admits a particular solution of the form 
3f(2 — tx) 

p 2 = 


(62) 


(63) 


(2 — 2 tx + t 2 x + t 2 x 2 ) 2 
We find that /i 2 also satishes the Erst equation given in (1571) and forms a compatible 
solution to the system of equations ([57]) . 

Substituting the expressions Ag and p 2 hr (1581) and integrating the resultant set of 
equations, 

3(2 tx 3 — t 2 x A + 2(1 + tx — t 2 x 2 )x — t 2 x 2 ) 


It = 


Ix = 


n = 


(2 — 2 tx + t 2 X + t 2 X 2 ) 2 
3(2 + t 2 x 2 — 4 tx — t 2 x ) 

(2 — 2 tx + t 2 x + t 2 x 2 ) 2 ’ 

3t(2 — tx) 


(2 — 2tx + t 2 x + t 2 x 2 ) 2 ’ ^ ^ 

we hnd 3J 2 = / 2 , where J 2 is given in flhOD . With the help of I\ and J 2 , we can derive 
the general solution for MEE. 


6. Hidden Symmetries 

An important application of Lie point symmetry is that it can be used to reduce the 
order of the underlying ordinary clifierential equation. It was observed that the order 
of the reduced ODE may admit more or lesser number of symmetries than that of the 
higher order equation. Such symmetries were termed as “hidden symmetries”. This 
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type of symmetry was first observed by Olver and later extensively investigated by 
Abraham-Shrauner and her collaborators [25, 26| EZ, 23 ESI 66]. 

The motivation for finding hidden symmetries of differential equations is the 
possibility of transforming a given ODE which has insufficient number of Lie point 
symmetries to be solved to another ODE that has enough Lie point symmetries such 
that it can be solved by integration. These hidden symmetries cannot be found through 
the Lie classical method for point symmetries of differential equations. 

A detailed study on the hidden symmetries of differential equations show that there 
can be two types of hidden symmetries. For example, if a n th order ODE is reduced in 
order by a symmetry group then two possibilities may occur. The reduced lower order 
ODE may not retain other symmetry group of the n th order ODE. Here the symmetries 
of the n th order equation is lost in the reduced equation. This lost symmetry is called a 
Type I hidden symmetry of the lower order ODE. Conversely, the lower order ODE may 
possess a symmetry group that is not shared by the n th order ODE. In this case, the 
lower order ODE has gained one symmetry. This is called as Type II hidden symmetry 
of the n th order ODE [25j E3 EH E3 EH [66] . 

Since we are focusing our attention on second-order ODEs, we again consider the 
MEE equation as an example and point out the hidden symmetries associated with this 
equation. 


6.1. Example: modified Emden equation 


Let us consider the MEE equation 
transformation 

y , 

X = -, t = z, 

y 


©. By introducing the following Riccati 


the MEE can be transformed to a linear third-order ODE, that is j-§ = 0. 


(65) 


The 


transformation (165]) is nothing but the invariants associated with the Lie point symmetry 
y-^ of the linear third-order ODE. In other words, the third-order equation = 0 has 
been order reduced to MEE by one of its point symmetry generator yj^. The other Lie 


point symmetries of the third-order linear ODE, = 0, are 


d 


d 


d 


d 3 y _ 

d 


Xl = *2 = 7T> X3 = ^ 7T> X4 = Z-, 

oz ay oy oz 


d 


d 


d 


Xb = z-z-, Xe = y-z~, X7 = 
oy oy 


yz-, 


z^d_ 

2 dz ' dy 

Substituting the transformation (1651) in the remaining vector fields given in 
they can be transformed into the following forms, namely 


( 66 ) 


* = I = 


Vo = 


-xe 


ii 


X 7T = V2 

ox 


Vs, 


d_ 
] dx 


J xdt ^ 

v 3 = 

(- 

t 2 \ . 

- xe 

dx ’ 


\x 

2 J 

, V 5 = ( 

--t] 

| xe~ 

- J xdt ^ 

’ 5 \ 

.x ) 


dx ’ 

Vj - Vs, 





— J xdt 


d_ 

dx' 


(67) 
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where i — 1, 2, 3,4, 5, 7, are the symmetry generators of the MEE (see Eqs. ffTUl) and 
(l89|) 1. While three of the vector fields {V \, V 4 and Vf) retain their point-symmetry nature, 
the remaining three vector fields (V 3 , V 3 and V 5 ) turns out to be nonlocal vector fields. 
All these vector fields satisfy the invariance condition and turn out to be the vector 
fields of the MEE. The local vector fields V\{— Vi), Va{— V 2 — V5) and Vt(— V 7 — Vq) 
match with the earlier ones (see Eq. (fTOjl ) whereas the nonlocal (V 3 , V 3 and V 5 ) vector 
fields emerge as new ones. 

Now we pick up Type-I and Type-11 hidden symmetries from them. As we pointed 
out in the beginning of this section, Type-II hidden symmetries of third-order ODEs 
are nothing but the symmetries gained by the second-order ODEs. The MEE admits 
eight Lie point symmetries (see Secj2]). In the above, we obtained only three Lie point 
symmetries of Eq.(Jl]). The remaining five Lie point symmetries are Type II hidden 
symmetries of the third-order ODE. These five symmetries can be gained from either 
non-local symmetries or contact symmetries of the third-order ODE. 

Type-I hidden symmetries of MEE are the symmetries which may not retain the 
symmetry group of the third-order ODE. In the present case, they turned out to be % 3 
and X 5 since these two vector fields cannot be found in (1671) . 

7. Nonlocal symmetries 

The study of hidden symmetries of ODEs brought out a new result. Besides point and 
contact symmetries, the ODEs do admit nonlocal symmetries (the symmetry is nonlocal 
if the coefficient functions £ and 7 ] depend upon an integral). The associated vector field 
is of the form V = f(t,x,J u(t, x)dt)-^ + r](t,x,J u(t, x)dt)j^.. Subsequently attempts 
have been made to determine nonlocal symmetries of ODEs. However, due to the 
presence of nonlocal terms, these nonlocal symmetries cannot be determined completely 
in an algorithmic way as in the case of Lie point symmetries. The determination of 
nonlocal symmetries for second-order ODEs was initiated by Govinder and Leach |28j . 
Their approach was confined to the determination of these nonlocal symmetries that 
reduce to point symmetries under reduction of order by t|. Later several authors have 
studied nonlocal symmetries of nonlinear ODEs [251, [26l [2711281 291 130]. Nucci and Leach 
have introduced a way to find the nonlocal symmetries [29J. In the following, we present 
a couple of methods which determine nonlocal symmetries associated with the given 
equation. We again consider MEE as an example in both the methods and derive the 
nonlocal symmetries of it. We also discuss the connection between A-symmetries and 
nonlocal symmetries. 

7.1. Method of Bluman et al. IWj 

In this method one essentially introduces an auxiliary “covering” system with auxiliary 
dependent variables. A Lie symmetry of the auxiliary system, acting on the space 
of independent and dependent variables of the given ODE as well as the auxiliary 
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variables, yields a nonlocal symmetry of the given ODE if it does not project to a point 
symmetry acting in its space of the independent and dependent variables. This method 
was first initiated by Bluman [32] and later extensively investigated by Gandarias and 
her collaborators [33| 34] [35]. 

Let the given second-order nonlinear ODE be of the form ([2]). To derive nonlocal 
symmetries of this equation, the authors introduced an auxiliary nonlocal variable y 
with the following auxiliary system [33, Ml [35], 

X - <j>(t,x,x) = 0, y = f{t,x,y). (68) 

Any Lie group of point transformation V = f(t,x,y)j^ + r](t,x,y + if(t,x,y)-j ^, 
admitted by (1681) yields a nonlocal symmetry of the given ODE (|2]) if the infinitesimals 
£ or rj depend explicitly on the new variable y, that is if the following condition is 
satisfied f 2 + rj 2 ^ 0. As the local symmetries of (168]) are nonlocal symmetries of (J2]) 
this method provides an algorithm to derive a class of nonlocal symmetries for the given 
equation. These nonlocal symmetries can be profitably utilized to derive the general 
solution for the given equation. Using this procedure, Gandarias and her collaborators 
have constructed nonlocal symmetries for a class of equations [33] [33., 35] • 

In the following, using the ideas given above, we derive nonlocal symmetries for the 
MEE. 

7.1.1. Example: modified Emden equation We introduce a nonlocal variable y and 
rewrite Eq. (JT]) in the form [35] 

x + 3xx + x 3 = 0, y — f(t,x,y ), (69) 

where f(t,x,y ) is an arbitrary function to be determined. Any Lie group of point 
transformation V — £(f, x, y)-j^ + rj(t,x,y)-j^ + if{t,x,y)-j^ admitted by (1691) yields a 
nonlocal symmetry of the ODE (ill) , if the infinitesimals £ and r) satisfy the equation 
e y +V 2 y^0. 

The invariance of the system (1691) under a one parameter Lie group of point 
transformations leads to the following set of determining equations £, r) and /, namely 

fn — 0, l[>x f £r — 6, T]xx fxiy 2 ‘If^xy T — 0, 

A + ffy - fit - fiy ~ M ~ fxV = 0, 

2 x 3 it + 2 fx 3 iy - y x x 3 + 3yx 2 + 3y t x + 3 fy y x + ryt + f 2 r] yy + ZfVyt + ftVy = 0, 

3 xit - itt ~ f 2 iyy - 2/ i yt + 3fxiy - ftiy + 3x 3 G + f x rjy + 2y tx + 2 fy xy + 3y = 0. (70) 

Solving the overdetermined system (1701) we obtain the following infinitesimal 
symmetry generator for the Eq. (]69]) : 

V = c(t)e y (x ^ + Jy), (71) 

with 

f(t,x) — X 


(72) 
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where c(t) is an arbitrary function of t. We note here that (T7T1) is not the only solution 
set for the determining equation (ITOj) . 

Solving the characteristic equation, we find two functionally independent invariants 
which are of the form 


x 


z = t, ( = —b x. 

X 


(73) 


In terms of these two variables, z and £, Eq.(pQ) reads as G + C 2 — 0. The general solution 

of this first-order ODE can be given readily as ( = - with k\ as an integration 

t + k i 

constant. Plugging this expression in the second equation in ( 1751 ) and rewriting it, we 
find 


This first-order ODE can be integrated straightforwardly to yield 

2 (t + k\) 
t 2 + 2 ki t — 2 k 2 ’ 

where k 2 is the second integration constant. Replacing k\ = I\ and —2 k 2 
we end up at the expression given in Eq. (lT5]h 


(74) 

(75) 

1 2 iii 


1.2. Connection between nojilocal symmetries and \-symmetries 


The exponential nonlocal symmetries admitted by (l69|) can also be derived from A- 
symmetries. To illustrate this we recall the following theorem from Ref. (4_0j. 


Theorem 1 Let us suppose that for a given second-order equation (d|) there exists 
some function f = f(t,x,x) such that the system [6E) admits a Lie point symmetry 
V = f(t,x,y)§- t + v(t,x,y )J^ + i>(t,x,y)-§^ satisfying Q + Vy + 0. We assume 
that z = z(t,x), ( = ((t,x,x) are two functionally independent functions that verify 
V(z) = 0, RT)(£) — o and are such that equation m can be written in terms of 

{z, (, C-} as a first-order ODE. Then 


1. The vector field V has to be of the form 

V = e Cy (f(t, x)d t + rj(t, x)d x + if(t, x, x)d y ) + C x d y , (76) 

where C and Ci are constants. 

2. The pair 


v = £,{t,x)d t + rj(t,x)d x , A = Cf. (77) 

defines a X-symmetry of the equation m and the set {z, (, f z } is a complete system 
of invariants of v PA 1 )]. 

With the choice C — 1, C\ — 0 and / = A, the vector field (1761) turns out to be 
V = e y (£(t, x)d t + r)(t , x)d x + x, x)d y ), 


( 78 ) 
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where £ and p are the infinitesimal coefficients of v and 0 = i fi(t,x,x ) satisfy the 
condition V ^ (0—A) | a = 0. This equation provides a linear first-order partial differential 
equation to determine 0 , that is 

0t + X A + 0*0 + 0A = A(£) A + £A 2 + v [a,( 1 ) 1 (A). (79) 

Let w = £<9i + rjd x be a A-symmetry of ([2j) for some A = A (t, x, x) and 0 = 0(t, x, x) 
be a particular solution of the equation (l79lh Then (1781) is a nonlocal symmetry of (121) 
associated to system (j68|) for / = A (t,x,cc) [40] • 

7.2.1. Example: modified Emden equation Using the above, we can demonstrate that 
the nonlocal symmetries found by Gandarias et.al for the MEE can be extracted from 
the A-symmetries themselves. To show this let us consider the A-symmetry /h -with 
A 3 = | — x (from Table 1). Substituting this expression in Eq. (1791) and solving the 
resultant partial differential equation we can obtain an explicit expression for 0 (t, x, x). 
Let us choose the simplest case ijj(t,x,x) — 0. In this case the left hand side of Eq. (l79j) 
disappears and the right hand side of it also vanishes automatically since it is nothing 
but the A-symmetry determining equation in which A 3 is a solution. Substituting A 3 = / 
in the second expression given in ( 1681 ) and integrating it, we fold 

y = log x — J xdt. (80) 

Now substituting the expressions £ = 0,77 = 1 ,0 = 0 and the above expression of y in 
d78|) . we obtain a nonlocal symmetry 

n 4 = xe -f xdt ^-. (8i) 

ox 

One can unambiguously verify that the vector field (1ST]) also satisfies the determining 
equation and turns out to be a nonlocal symmetry of the MEE. We mention here that 
the nonlocal symmetry (IHTj) had already been observed in the order reduction procedure 
(see Eq. (I67]P . 

The other choices of A and/or fi>(t,x,x) will generate new nonlocal symmetries for 
the MEE. For example, the choice 0 = c(f),£ = 0, r/ = c{t)x and A = — x — —, provides 
another nonlocal symmetry (1711) through the above said procedure. In this way, one can 
also construct nonlocal symmetries from the A-symmetries. 

7.3. Method of Gladwin Pradeep et.al 

In a recent paper Gladwin Pradeep et.al proposed yet another procedure to determine 
the nonlocal symmetries for the given equation [31j. In the following, we briefly recall 
the essential ideas behind this method with reference to MEE. 

The MEE (JT() can be transformed to the second-order linear ODE = 0 through 
the nonlocal transformation u = xe-f xdt . To explore the nonlocal symmetries associated 
with (H]) , the authors used the identity „ — x + x (which can be directly deduced from 
the nonlocal transformation u = xe-f xdt ) [3lj. This nonlocal connection between the 
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free particle equation and MEE allows one to deduce the nonlocal symmetries of Eq. 
dU) in the following manner. 

Let £ and r) be the infinitesimal point transformations, that is u' = u + er](t,u), 
T = t + e£(£, u), associated with the linear ODE ^ = 0. The symmetry vector 
held associated with this infinitesimal transformations reads as V = £y + rand 
its first extension is given by Pr^V = £y + r]-^ + (?) — Then we denote the 

symmetry vector field and its first prolongation of the MEE ([T]) as fl = 5j^ + and 
Pr^Q = Sijj + [j,-^ + (/i — x8)-jjr, where <5 and /i are the infinitesimals associated with 
the variables t and x, respectively. 

From the identity - — - + x, the authors defined a new function, say X 

- = - + x = X. (82) 

u x 

In the new variable X, the MEE turns out to be the Riccati equation, that is 

X + X 2 = 0. (83) 


The symmetry vector field of this equation can be obtained by using the relation 


X — - and rewriting V 1 = £-4 + TJjj- + (?) 


Oik as 


i/I =4 + 




u 


u 


d_ 

JdX 


= E. 


(84) 


We note that Eq. (1831) . being a first-order ODE, admits infinite number of Lie point 
symmetries. These Lie point symmetries of Eq. (T83|) become contact symmetries of the 
linear second-order ODE ^ = 0 through the relation X — 

Similarly one can rewrite O 1 = + /iy: + (ji-x 8 )Jp-, using the relation X = - + x, 


as 

nl =4 + (<■?* +i)m+(a ■ *4) ix= s - (85) 

As the symmetry vector fields £ and 5 correspond to the same equation (155|) . their 
infinitesimal symmetries must be equal. Therefore, comparing equations (1841) and (1851) 
one obtains 


£ = h, -x£ = {— \x + l)/i + ji~. (86) 

U U z X z X 

Rewriting the second equation given in (I86p we can obtain the following first-order ODE 
for the unknown function /a, that is 


-fi + (~\x + l)fi = - x£. (87) 

x x z at u 

The free particle equation (yr = 0 admits eight Lie point symmetries which are 
given in Eq. flTZjh Substituting these symmetries (£j, iji), i = 1,2,, 8, and u = xef xdt , 
in Eq. (1871) . we get the following seven first-order ODEs for /a, 

— a + (1 —— x)ia + (x 2 + x)x~ 2 e~f xdt = 0, (88a) 

x x l 

—\a + (1--x)/i — (x — tx 2 — tx)x~ 2 e~ $ xdt = 0, (886) 

X X 1 
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-H + (1 - 7^x)[i = 0 , 

x x- 

11 r 

— ii + (1--x)n + (x 2 + x)xe~ J xdt = 0, 

X x^ 

1 . 1 . 

~(l + (1-ttx)/x + X = 0, 

X x^ 

1 1 

—jx + (1--x)/i + 2 tx — 1 = 0, 

X X" 

—/i + (1 —4-x)/x + 2(x 2 + x)x 2 e 2 f xdt — 1 = 0. 
x x l 


(88 c) 
(88d) 
(88 e) 
( 88 /) 
( 880 ) 


Integrating each one of the above first-order linear ODEs we can obtain the function ji. 
Substituting the symmetries S(= £) and ^ in 12 = 5-^ + we get following nonlocal 
symmetries, 

d 


Ql ~ dt’ Q ’ 2 ~ (x \ 


fl.-i — 


t 

X 


d 


xe ~f xdt — 

ox 
d 


~f xd t Q = xe -f xdt 


d 


dx ’ 4 dx ’ 

O, = xef xdt Z - ( f x(x + x 2 )el {2x)dt d?j xe~f xdt 


dt 


ft 6 = tJL- xe-f (xfx)dt ^ xe-( xdt dtj 


Q r = t 2 ^ + xe~f xdt (^j (1 - 2 tx)e-f xdt dt^ 

= txef xdt + xe~$ xdt (f ( tx 3 + (tx - l)x)e/ (2x)dt df 


_d_ 

dx 


(89) 


of equation ©■ One may observe that some of the nonlocal vector fields 02,03 
and 0 4 had already been found as hidden symmetries earlier. It is a straightforward 
exercise to check that all these nonlocal symmetries satisfy the invariance condition 
5^ + — ytb 2 ) = 0 , where and /b 2 ^ are the first and second prolongations. 

We mention here that these nonlocal symmetries can also be related to A-symmetries 
through the theorem given in Sec. 17.21 

To derive the general solution of the given nonlinear ODE one has to solve the 
Lagrange’s system associated with the nonlocal symmetry. For the vector field 0 4 , the 
underlying equation reads (Eq. (j89jl ). 

dt dx dx . . 

—— = — = —-—. (90) 

0 x x — x 2. 

Integrating Eq. ([90]h we find u — t and v — - + x. Following the procedure described in 
Sec l4.31 we can obtain the general solution of (JT]) as in the form Eq. (fT5|) . 


8. Telescopic vector fields 


Telescopic vector fields are more general vector fields than the ones discussed so far. 
The Lie point symmetries, contact symmetries and A-symmetries are all sub-cases of 
telescopic vector fields. A telescopic vector field can be considered as a A-prolongation 
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where the two first infinitesimals can depend on the first derivative of the dependent 
variable [32 , 30] , In the following, we briefly discuss the method of finding telescopic 
vector fields for a second-order ODE. We then present the telescopic vector fields for 
the MEE. 

Let us consider the second-order equation (EJ) . The vector field 


„(*)=£» + » +{ (D» +{ W» 


’ dt Ox 


Ox 


dx 


is telescopic if and only if 


(91) 


Z=Z(t,x,x), i) = t l(t,x,x), C (1) = C (1> (t,x,±) 


(92) 


with £( 2 ) is given by 

C ( 2 ) = D[cW] _ 0jD[e] + C (1) + X D[£\ — D[r]\ ( ^.(i) _ (93) 

n - 

To prove that the telescopic vector fields are the more general vector fields, let us 
introduce two functions g\ and g 2 in the following form, namely 

= Ja(t,i,i) = (94) 

g — xt, g — x£ 

We can rewrite the prolongations (A 1 ) and £( 2 ) using the above functions gi and g 2 as 
follows: 


c (1) = D[rj\- xD[£\ + (g x + g 2 (fi){g - x£), (95) 

C (2) = D[ C (1) ] - 0xD[£] + ( 9l + g 2 m (1) - <K)- (96) 

The relationship between telescopic vector fields and previously considered vector fields 
can be given by the following expressions [32, 130] 

C (1) =77 (1) + (9i + 924 >)(.V-b£), (97) 

c (2) = V {2) + (.01 + 32</>XC (1) - <K)- (98) 


In the above vector fields if we choose g\ = g 2 = 0 and + g 2 = 0 we get the Lie 
point symmetries. The choice g\ = g 2 = 0 and £| + g 2 X 0 gives the contact symmetries. 
To get A-symmetries, we should choose g\ X 0 and £| + r/J = 0. As a consequence it can 
be considered as the more general vector field. 


8.1. Example: modified Emden equation 

To End the telescopic vector fields admitted by the MEE equation, we have to solve the 
invariance condition v^ 2 \fi) = + + + C^§f = X with £ and g are functions 

of (t, x, x) and (3 1 ) and are defined through (j97j) and (BHD respectively. Substituting 
Eq. (]TJ) in the invariance condition, we obtain 

- (3x + 3x 2 )g - 3xC (1) - C (2) = 0. (99) 

Solving equation (l99lh we obtain a telescopic vector field which is of the form 


( x ') 

1 d +1 

/ x 2 — X \ 

A +1 

/ 6xx A 


V (x 2 + x) 2 ) 

ox 

V (x 2 + x) 2 / 

I • I 

Ox 

V(x 2 + x) 2 / 

Ox 


7i = 


( 100 ) 
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6 = 0 , r )i = 


X (i) _ X 2 - X ( 2 ) _ Qxx 

( 2 i •\2 5 / 2 i •\2 5 / o i • \2‘ 

+ X) + X) ( X z + X) 


A second telescopic vector held is found to be 


72 = - 


t( 2 — hr) 




5 


+ 


f 2 (x — x 2 ) + 4te - 2 \ 9 
(f 2 (x 2 + x) — 2 to + 2) 2 y cte V (t 2 (x 2 + x) — 2 tx + 2) 2 ) dx 

6 (tx — l)(ti + x) \ d 
(; t 2 (x 2 + i) — 2 te + 2) 2 / dx 


and its components are given by 

t( 2 — tx) 


6 = 0 , 772 = 


d 2) = 


(t 2 ( x 2 + x) ■— 2tx + 2 ) 
6(te — l)(hr + x) 


A (i) t 2 (x — x 2 ) + 4te 

2 ’ 82 = 


(t 2 (x 2 + x) — 2tx + 2 ) 


(t 2 (x 2 + x) — 2 tx + 2 y 


( 101 ) 


( 102 ) 


2 1 


(103) 


The invariants associated with a telescopic symmetry vector held can be derived 
by solving the associated characteristic equation. For the vector held 71 , it reads 
dt dx dx 

O’ “ “ ~ 


(104) 


(x 2 +x) 


(x 2 -\-x) 


Using the procedure discussed in Sec l2.2,ll we can integrate the above characteristic 
equation to obtain the integral given in Eq. (j59li . Repeating the procedure for the second 
telescopic vector held (I102[) we end up at the second integral given in Eq. (j60]) . From 
these two integrals we can derive the general solution of (ITT) . 


9. Conclusion 

In this paper, we have reviewed continuous symmetries of second-order ODEs and 
elaborated the methods of finding them. To begin with we have considered Lie point 
symmetries and presented Lie’s invariance analysis for a second-order ODE. To illustrate 
the method, we have considered the modified Emden equation (MEE) as an example. We 
have also discussed few applications of Lie point symmetries. We have demonstrated the 
connection between symmetries and conservation laws by recalling Noether’s theorem. 
Few conserved quantities including energy have been identihed for the MEE through this 
theorem. We then considered the velocity dependent transformations and presented the 
method of finding contact symmetries for the second-order ODEs. We have also pointed 
out the contact symmetries of the MEE. We have also recalled hidden symmetries 
of the MEE. Some of them are found to be exponential nonlocal symmetries. The 
connection between symmetries and the integrating factors of ODEs was discussed 
through A-symmetries approach and adjoint symmetries method. The method of finding 
A-symmetries, adjoint symmetries, integrating factors and their associated integrals of 
a second-order ODE are discussed elaborately and illustrated with MEE as an example. 
We have also pointed out the connection between exponential nonlocal symmetries and 



















Symmetries 


29 


A-symmetries. Finally, we have considered a more generalized vector field, namely 
telescopic vector field and discussed the method of finding these generalized vector 
fields. For the MEE we have also brought out a couple of telescopic vector fields. We 
have also derived the general solution of MEE from each one of these symmetries. The 
symmetry methods presented here are all extendable to higher order ODEs. Through 
this review, we have emphasized the utility of symmetry analysis in solving ODEs. 
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